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Abstract —In this paper, we propose an end-to-end channel 
model for molecular communication systems with metal-oxide 
sensors. In particular, we focus on the recently developed table 
top molecular communication platform. The system is separated 
into two parts: the propagation and the sensor detection. There is 
derived, based on this, a more realistic end-to-end channel model. 
However, since some of the coefficients in the derived models are 
unknown, we collect a great deal of experimental data to estimate 
these coefficients and evaluate how they change with respect to 
the different system parameters. Finally, a noise model is derived 
for the system to complete an end-to-end system model for the 
tabletop platform. 

Index Terms —Molecular communication, channel model, 
tabletop platform, metal-oxide sensor. 

I. Introduction 

In molecular communication, information is transferred via 
chemical signals. Unlike the traditional wireless communica¬ 
tion systems, molecules are used to physically carry infor¬ 
mation [1], The transmitter releases molecules, often called 
messenger/information molecules, that propagate through dif¬ 
fusion [2]-[4], medium flow [5], [6], or active transport [7]- 
[9] to arrive at the receiver. The information can be encoded 
in several ways. For example, different information can be 
represented by different concentration, number, and/or type of 
molecules [6], [10], [11]. 

In the literature, most prior work has focused on the 
communication via diffusion, or diffusion assisted by flow. 
Most of these works, moreover, are only theoretical-based 
analyses of these systems. The authors in [12] analyzed the 
maximum data transmission rates of diffusion-based systems. 
There has also been several work regarding proper modu¬ 
lation techniques for molecular communication [12]. Some 
potential applications have been addressed in [13], and one 
of the promising applications for molecular communication is 
considered to be the biomedical field. 

The first experimental platform for molecular communica¬ 
tion was presented in [14]. This system, which was designed to 
be inexpensive and flexible, proved the feasibility of molecular 
communication, especially on the macroscale. This testbed has 
an electronically controllable spray that acts as a transmitter, 
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a metal-oxide gas sensor that acts as a receiver, and a fan to 
assist the propagation. Two Arduino microcontrollers are also 
used to control the actions of the transmitter and the receiver. 
Fig. 1 highlights this system and its different components. 

One of the main advantages of molecular communication is 
its multi-scale property: molecular communication can be used 
on both the microscale and macroscale. Although the tabletop 
platform designed in [14] is a macroscale system, some of 
the components can be shrunk to micro- and nanoscale. 
For example, metal-oxide sensors can be easily shrunk to a 
nanoscale, and be used for detection of different biological 
compounds [15], [16]. Therefore, understanding how these 
sensors work can be very beneficial for both macroscale and 
microscale molecular communication. 

As it was shown in [17], the system response of the testbed 
departs from the previous theoretical channel models used 
in the literature. Moreover, it was shown that the system 
tends to be nonlinear, where the nonlinearity can be modeled 
as noise. In [17], we attempted to add correction terms to 
the theoretical channel models. It was demonstrated that the 
proposed corrections to the system response would agree with 
experimental results. That work, however, considered only 
a single separation distance between the transmitter and the 
receiver and a single spray duration. 

In this work, we build on our previous work and derive a 
more realistic analytical model for the testbed. The following 
are the new contributions in this work (compared to [17]). 

• We consider a system model where the sensor and particle 
propagation are two separate systems. Previous work has 
modeled the tabletop platform as one system. By using 
two separate systems connected in series, a more realistic 
model is achieved. 

• We consider the end-to-end system response and how 
it changes with different system parameters such as 
distance, spray duration, and initial voltage (i.e., the 
initial concentration of the messenger molecules in the 
environment). Our previous work has considered a system 
model based only on a fixed set of system parameters. 

• A newly defined additive noise model is also presented 
for this systems, and is analyzed. 

The rest of the paper is organized as follows. Section II 
defines the system models, and Section III analyzes the effect 
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Fig. 1: The tabletop platform with the MQ-3 metal-oxide 
sensor. 
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Fig. 2: Block diagram of an end-to-end system. 


of the system parameters on channel response. A universal 
channel model is established in Section IV, and Section V 
concludes the paper. 


II. System Model 


As shown in Fig. 1 the transmitter is a spray, and the receiver 
is a metal-oxide sensor. A fan is used to assist the propagation, 
and the isopropyl alcohol is used to carry information. The 
end-to-end system block diagram of the tabletop platform is 
presented in Fig. 2. Essentially the system input is a spray 
pattern, and the system output is the voltage readings. The 
system has two components: the propagation module, where 
the sprayed alcohol is carried to the sensor, and the sensor 
module, where the detected alcohol concentration is converted 
to voltage. 

An impulse input is generated through a very short and 
instantaneous spray (e.g., 100ms). The system response to the 
impulse input is the impulse response of the system. In our 
previous work [17], the end-to-end system response of the 
system was based only on the chemical propagation model. 
Specifically, the model considered was 


C(t) = -j= exp ( -b 


(d — ct ) 5 


(1) 


where a, b, and c are correction factors introduced to match 
the theoretical and experimental results, and d is the separation 
distance between the transmitter and the receiver. 

In this work, we consider a more accurate model by 
incorporating the sensor block as well. 
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Fig. 3: The sensor anatomy. 


gasses as well as volatile liquids such as alcohol [18], They 
can also be shrunk to nanoscale and be used to detect vari¬ 
ous biological compounds [15], [16]. Moreover, metal-oxide 
sensors can be used to detect different chemicals [19]. 

The MQ-3 is a thin film tin dioxide (SnC> 2 ) sensor, where an 
aluminium oxide (AI 2 O 3 ) tube is covered with a thin layer of 
SnC> 2 . There is also a heater coil inside the aluminium oxide 
made from nickelchromium alloy, which is used for resistive 
heating. Fig. 3 shows the anatomy of the MQ-3 sensor. 

The sensor works as follows. First, the heater coil heats 
up the sensing layer (i.e., the SnC >2 layer). According to the 
sensor’s user manual, to achieve the optimal sensor sensitivity, 
the sensor must be heated 24 to 48 hours prior to use. When 
the SnC >2 layer heats up, it becomes a semiconductor. When 
alcohol vapour approaches the sensor, it will go through an 
oxidization reaction that will in turn change the resistance of 
the SnC >2 sensing layer [20]. In particular, the resistance of the 
SnC >2 sensing layer decreases as the concentration of alcohol 
increases in the vicinity of the sensor. 

There is a well-known relationship between the resistance 
of a metal-oxide sensor and the concentration of target 
molecules [ 21 ] 

R s “ ai C n (2) 

where R s is the resistance of the sensor; C is the concentration 
of molecules, and a\ and n are constants. In (2), the constant, 
a±, is empirically chosen, and the power law exponent, n, is 
specific to the type of target molecule. It can also be calculated 
theoretically from the reduced depletion depth and reduced 
reactivity [ 21 ]. 

To estimate the constant n provided in (2) for the MQ-3 
sensor, we use the sensitivity graph in the MQ-3 data sheet. 
In this graph, shown in Fig. 4, the ratio R s /Rq is plotted 
against the concentration of different target gases, where Rq 
is the resistance of the sensor at a specific concentration of the 
target gas Cq. In this plot, both axes are in log scale. Therefore 
we have 

log - lo S (o- 2 C{t) n ), (3) 

where 02 = ai /Ro is a new constant. From (3) it is clear that 
n is the slope of the line corresponding to alcohol in Fig 4. 
Therefore, we estimate the value of n = —0.65, and the sensor 
resistance is related to the concentration with 

R s (t) = ai C(t)-° S5 , (4) 

where a \ is an unknown constant. 


A. Metal-Oxide Sensors 

The sensor used in the platform is the MQ-3 metal-oxide 
sensor. Metal-oxide sensors are cheap and can detect various 


B. End-to-End Response 

As shown in Fig. 2, the system consists of propagation 
and sensor detection. For the propagation phase, the expected 
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Fig. 4: The sensor sensitivity plot from the datasheet [21], 


TABLE I: System Parameters. 


Distance (m) 

2 

3 

4 

5 

Spraying Duration (ms) 

50 

100 

150 

200 

Initial Voltage (V) 

1.0 

1.3 

1.6 

1.9 


system response is given by [22]: 


C{t) = M 


V 47t Dt 3 


exp 


( ( d-vt) 2 \ 

V 4£>t ) ’ 


(5) 


where M is the number of alcohol molecules released by the 
short spray; d is the distance between the transmitter and the 
receiver; D is the effective diffusion coefficient, and v is the 
average velocity of the wind flow. Substituting (5) into (4), 
the end-to-end system response is given by 


R s {t) = ai 


M — exp 
\J AnDt 3 


0 d-vt) 2 \ 
ADt ) 


( 6 ) 


The change in sensor resistance is measured using a simple 
voltage divider circuit. 

Since M, a±, the effective diffusion coefficient D, and the 
average velocity v cannot be measured accurately, the end-to- 
end system response of the system is given by 


h{t\ a , b,c) = a 


d ( 

7^ exp y 


(d — ct) 2 


-0.65 


(7) 


where a, b, c are unknown coefficients, and d is the distance 
between the transmitter and the receiver. Specifically, a is 
a the product of constant m and the number of transmitted 
molecules M; b is the effective diffusion, and c is the average 
velocity. By obtaining the proper values for the coefficients, 
a generalized universal model can be established. In the next 
section, we use experimental data to estimate the values of 
these three coefficients, and observe how they change with 
respect to different system parameters. 


III. Universal Model Based on Experiments 

To find the value of the unknown coefficients as a function 
of system parameters, the tabletop platform shown in Fig. 1 
is used to measure the end-to-end impulse response of the 
system using, as messenger molecules, isopropyl alcohol. To 
see the effects of different system parameters on the channel 
responses, three system parameters are considered: distance 
between the transmitter and the receiver, spraying duration 
(which is related to the number of the transmitted molecules), 
and initial voltage or initial concentration of the messenger 
molecules in the environment. Each parameter is varied over 
four different values as shown in Table I. Therefore, there 
are 4x4x4=64 different cases, and the end-to-end impulse 
response of each case is measured during 10 different trials. 
Thus, to collect all the data to estimate the unknown co¬ 
efficients, in this paper, we carry out, in total, 640 sets of 
experimental trials. 

A. Coefficients 

To estimate the value of the unknown coefficients in the 
analytical end-to-end model (7), we use a nonlinear least 
squares curve fitting technique used in our previous work 
[17], Assuming that there are N points in the sensor mea¬ 
surement from each trial, then (7) can be discretized as 
h(tk, p) (k £ {1,2, ••• ,N}), where tk are sampling times, 
and p = [a, 6, c] T is the parameter vector representing the 
three unknown constants. The coefficient estimation problem 
can be formulated as 

N 

arg min E (0(t k )-h(t k , p)) (8) 

p fc =i 

where 0{tk ) are the experimental observations. 

The coefficients are estimated for all the 640 data sets 
using the least squares estimation technique. For example, 
the fitted results are shown in Fig. 5 for the case when the 
distance between the transmitter and the receiver is 2 m, the 
spray duration is 150 ms, and the initial voltage is 1.3 volts. 
It is apparent that our model function fits very well to the 
measurement data. The root mean square error (RMSE) is 
0.0484, which is a bit smaller than with the previous model. 

Fig. 6 shows the effect of the system parameters on the 
coefficients. Each point in the plot represents the average 
value of the corresponding coefficient across 160 trials. The 
standard deviation bars represent the standard deviation across 
trials. As can be seen, coefficient a changes frequently with 
changing system parameters, while coefficient c varies slightly 
with respect to system parameters, and b is mostly a con¬ 
stant. Specifically, coefficient a changes with distance, spray 
duration, and initial voltage; c changes only with distance and 
initial voltage (not with spraying duration). These results agree 
with the physical characteristics of the system. For example, 
the effective diffusion coefficient is expected to be a constant. 
From Fig. 6(b), we can see that as the spray duration increases, 
a increases linearly and the variance of a also increases. This 
is because more molecules are released as the spray duration 
increases thereby increasing M in (6). 
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Fig. 5: One fitting example with three unknown coefficients when system parameters are 2 m, 150 ms, and 1.3 V. 





(a) Distance effect on the coefficients. (b) Spraying Duration effect on the coefficients. (c) Initial Voltage effect on the coefficients. 

Fig. 6: Averaged effect of the system parameters. 
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(a) The averaged effect of system parameters on (b) The averaged effect of system parameters on (c) The averaged effect of system parameters on 
coefficient a. coefficient c. the standard deviation of the additive noise, a. 

Fig. 7: The averaged effect of system parameters to the coefficients, a and c, and standard deviation of the additive noise, er. 
The x-axis indicates the index of the system parameter vectors as represented in the legend. 


Based on these observations, we assume b to be constant that since the variance of b is so small, this is a very good 

and be equal to the average values of b over 640 trials. Note estimation for this coefficient. Let b* be this average value 
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with 6* = 0.1950. 

When we perform least squares curve fitting with the fixed 
coefficient b* and unknown coefficients a and c, the fitted 
curves are as good as the ones shown in Fig. 5. Therefore, we 
leave out those plots to avoid duplications. Fig. 7(a) and 7(b) 
show the behaviour of the coefficients with the fixed b*. As 
can be seen, both a and c change with respect to the system 
parameters. Let d, s, and is be the distance, spraying duration, 
and the initial voltage, respectively. Then from (7) we have 


hit - aM - C)= "[VMr¥ 
a = f(d,s,v), c = g(d,is), 
h(t ; d, s, v) = 


exp 


( _ (d - ctf \ 
V 4 b*t ) 


-. -0.65 


(9) 


f(d, s, is) 


yj An r6*i 3 


exp ' 


(d-g(d, is)tf 

46* i 


-0.65 


where h(t ; d, s , is) is a function of time, t, with three coeffi¬ 
cients. 


Since f(d, s, is) is obtained from the average value of coeffi¬ 
cient a over different trials, the mean of N becomes zero for 
all cases. As represented in Fig. 7(c), however, the standard 
deviation of N changes with system parameters showing 
behaviour similar to a. Thus, it is also possible to express 
the standard deviation as a function of the system parameters 
as in analyzing coefficient a. 

E\N] = 0 

1—1 - (14) 

\JVar pV] = a = L(d , s, is) 

where E\-\ and Var\-\ denotes the expectation value and 
variance, respectively. 

IV. A Complete End-to-End Model 

In this section we derive a complete end-to-end model 
for the tabletop platform. From Section III, coefficient a, its 
additive noise’s variance, and c have a close relationship with 
the system parameters. Therefore, our goal is to estimate the 
functions f(d , s, is), g{d, is), and L(d, s, is). 


B. Noise 

An important part of the system is the random nature of 
the end-to-end impulse response. Some of the factors that 
contribute to the randomness are: the spray, which is not 
precise enough to spray the same amount of alcohol across 
trials, the random propagation due to diffusion and turbulent 
flows, and other phenomena. Therefore, to obtain a complete 
system model, the random effect must be represented as noise. 
To establish a noise model, an additive noise is defined as 
the difference between the experimental observation and the 
model function shown below. 

Ni(t) = Oi(t) - h(t;d,s,is) (10) 

where N, t (t) represents the additive noise and 0,(f) represents 
the observation or measurement for each trial i (i = 1, 2, 10 
for a specific distance, spray duration, and initial voltage). 
From Fig. 7(a) and 7(b), we can see that for coefficient c, the 
standard deviation across the trials is slight. Therefore, since 
only coefficient a has a large variance across different trials 
we have 


Oi(t) = 

(f(d,s,is) + N) 

where N is the noise introduced by the system, with 


d ( 

<N 

50 

IT 

I 

s/Antb*t 3 P V 

46* i 7 J 


-0.65 


( 11 ) 


Ni(t) = N 


VAnb*t 3 


exp - 


(d — g(d, is)t f 


46* t 


-0.65 


( 12 ) 


Let aoi be the best least squares fit for coefficient a for the 
trial experimental observation Oi(t). Then a noise sample for 
N is obtained by: 


N = a 0i - f(d,s,is) (13) 


A. Estimating the Coefficient Functions 

From Figs. 7(a), 7(b), and 7(c), it can be seen that the coef¬ 
ficients a, c and the noise variance change relatively linearly 
with respect to the different system parameters. Therefore we 
estimate the functions f(d, s, is), g{d, is), and L{d, s, is) as 
linear functions, respectively 

f{d,s,is) = /3 { d t) d +/3 { s f) s +/3l f) is +f3 ( 0 f) , (15) 

g{d,u)=^d + ^v + ^\ (16) 

L(d, s, is) = f3 { d L) d + pi L h + tfpv + /3q L \ (17) 

where all the (3 variables are constants. To estimate the values 
of the /3 constants, we use again the linear least squares 
method. Table II shows the estimated values of these constants. 

TABLE II: The estimated /3 coefficients in (15)—(17). 


Coefficient 

Value 

Coefficient 

Value 

Coefficient 

Value 


-0.4188 

A 3) 

0.0709 

ti L) 

-0.1258 

d f> 

0.0098 


- 

d L) 

0.0014 


-1.7873 

dl 9) 

0.1362 

pl L) 

-0.2403 


4.6469 


-0.0427 


0.9738 


B. Verification 

To verify the estimation in the previous section, two data 
sets are selected: 1) a distance of 2.5 m, a spray duration 
of 130 ms, and an initial voltage of 1.5 V and 2) 3.5 m, 
170 ms, and 1.1 V. Fig. 8 compares the comparison between 
the estimated system model and the averaged observations, 
and it shows that the two results match fairly well. Thus, the 
newly developed end-to-end channel model can be utilized for 
whichever values of the distance, spraying duration, and initial 
voltage. 
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Fig. 8: Comparison of the estimation from the obtained model 
and the experimental channel observations. 

V. Conclusions 

In this work, we developed an end-to-end channel model 
based on experimental observations using the tabletop plat¬ 
form. Unlike previous work on the channel model, it is repre¬ 
sented as two separate systems of the sensor and particle prop¬ 
agation making it a more realistic model. Also, it is generalized 
with three different system parameters: distance, spraying 
duration, and initial voltage. Furthermore, the randomness of 
the system is represented as an additive noise model. The 
two data sets also confirm that the newly developed universal 
channel model works well. For future work, we will propose 
new modulation techniques appropriate for this platform using 
this channel model, and compare the performances depending 
on the system parameters. 
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